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Abstract 

We develop a universal approach to the one-loop effective field theory (EFT) using the Co¬ 
variant Derivative Expansion (CDE) method. We generalise previous results to include broader 
classes of UV models, showing how expressions previously obtained assuming degenerate heavy- 
particle masses can be extended to non-degenerate cases. We apply our method to the general 
MSSM with non-degenerate stop squarks, illustrating our approach with calculations of the 
coefficients of dinrension -6 operators contributing to the hgg and /177 couplings, and comparing 
with exact calculations of one-loop Feynman diagrams. We then use present and projected 
future sensitivities to these operator coefficients to obtain present and possible future indirect 
constraints on stop masses. The current sensitivity is already comparable to that of direct LHC 
searches, and future FCC-ee measurements could be sensitive to stop masses above a TeV. The 
universality of our one-loop EFT approach facilitates extending these constraints to a broader 
class of UV models. 
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1 Introduction 


In view of the overall consistency between the current measurements of particle properties and 
predictions in the Standard Model (SM), a common approach to the analysis of present and 
prospective future data is to describe them via an effective field theory (EFT) in which the 
renormalizable SM d = 4 Lagrangian is supplemented with higher-dimensional terms composed 
from SM fields |1,2|. To the extent that this new physics has a mass scale that is substantially 
higher than the energy scale of the available measurements [3j, the EFT approach is a powerful 
way to constrain possible new physics beyond the SM (BSM) that is model-independent 4]|6]. 
The d = 6 operators in this Effective SM (ESM) were first classified in |l[^J with a complete 
basis using equations of motion to eliminate redundancies |2j being first presented in |9j. There 
have been many studies of various aspects of these dimension-6 operators [^J and a short review 
can be found in 
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The EFT approach may well be a good approximation if the new physics affects precision 
observables at the tree level, or if it is strongly-interacting. In these cases the new physics mass 
scale is likely to be relatively high, and considering the lowest-dimensional EFT operators may 
well be sufficient. However, the EFT approach may have limitations if the new physics has 
effects only at the loop level, or is weakly interacting. In these cases, the EFT approach may 
be sensitive only to new physics at some relatively low mass scale, and the new physics effects 
may not be characterised well by considering simply the lowest-dimensional EFT operators. 

Examples in the first, ‘safer’ category may include certain models with extended Higgs 
such as two-Higgs-doublet models, or some composite models. Examples in the 
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sectors 

second category may include the loop effects of supersymmetric models. However, even in this 
case it is possible that precision electroweak and Higgs data may provide interesting constraints 
on the possible masses of stop squarks, which have relatively large Yukawa couplings to the 
SM Higgs field. In particular, the EFT approach may be useful in the framework of ‘natural’ 
supersymmetric models with stops that have masses above 100 GeV but still relatively light 
compared to other supersymmetric particles. 

Important steps towards the calculation of loop effects and the simplification of their match¬ 
ing with EFT coefficients have been taken recently by Henning, Lu and Murayama (HLM) [12, 
13 . In particular, they use a covariant-derivative expansion (CDE) [17,18 to characterise new- 
physics effects via the evaluation of the one-loop effective action. They apply these techniques 
to derive universal results and also study some explicit models including electroweak triplet 
scalars, an extra electroweak scalar doublet, and light stops within the minimal supersymmet¬ 
ric extension of the SM (MSSM), as well as some other models. They also discuss electroweak 
precision observables, triple-gauge couplings and Higgs decay widths and production cross sec¬ 
tions 113] , and have used their results to derive indicative constraints on the basis of present 
and future data |l2|. 

In this paper we discuss aspects of the applicability of the EFT approach to models with 
relatively light stops, exploring in more depth some issues arising from the work of HLM [12jl3|. 
As they discuss, using the CDE and the one-loop effective action is more elegant and less time- 
consuming than a complete one-loop Feynman diagram computation. On the other hand, they 
applied their approach to models with degenerate soft supersymmetry-breaking terms for the 


1 This EFT approach that we follow, in which the SU( 2)l x U(\)y electroweak symmetry is linearly 
realized, is to be distinguished from a non-linear EFT based on the chiral electroweak Lagrangian |7 
and the mor e ge neral anomalous coupling framework of a C/(T)em effective Lagrangian |8j. 


2 See 0§|lO for some examples of earlier work and [dj-jbj 9 
studies. 
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for a sampling of more recent 
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stop squarks, and we show how to extend their approach to the non-degenerate case, with specific 
applications to the dimension -6 operators that contribute to the hgg and /177 couplings. Our 
extension of the CDE approach would also permit applications to a wider class of ultra-violet 
(UV) extensions of the SM and other EFT operators. 

Another important aspect of our work is a comparison of the EFT results with the corre¬ 
sponding full one-loop Feynman diagram calculations also in the non-degenerate case, so as to 
assess the accuracy of the EFT approach for analysing present and future data. 

In a recent paper, together with Sanz, two of us (JE and TY) made a global fit to dimension -6 
EFT operator coefficients including electroweak precision data, LHC measurements of triple¬ 
gauge couplings, Higgs rates and production kinematics | 6 j. Here we use this global fit to 
constrain the stop mass nq- and the mixing parameter Xt, comparing results obtained using 
the EFT with those using the full one-loop diagrammatic calculation. The bounds on rrq 
and Xt are strongly correlated, and we find that the EFT approach may yield quite accurate 
constraints for the limits of larger mi and Xt. However, there are substantial differences from 
the full diagrammatic result for smaller m1 and Xf. In this case the diagrammatic approach 
gives indirect constraints on the stop squark that are quite competitive with direct experimental 
searches at the LHC. We also explore the possible accuracy of the EFT for possible future data 
sets, including those obtainable from the LHC and possible e + e - colliders [^] For example, 
possible FCC-ee measurements 22 may be sensitive indirectly to stop masses > 1 TeV. 

The layout of this paper is as follows. In Section 2 we introduce the covariant derivative 
expansion (CDE) and discuss its application to the one-loop effective action, highlighting how 
the HLM approach 12,13 may be extended to the case of non-degenerate squarks. As we 


discuss, one way to achieve this is to use the Baker-Campbell-Hausdorff (BCH) theorem to 
rearrange the one-loop effective action, and another is to introduce an auxiliary expansion 
variable. Results obtained by these two methods agree, and are also consistent with the full 
one-loop Feynman diagram result presented in Section 3. Analyses of the current data in 
the frameworks of the EFT and the diagrammatic approach are presented in Section 4, and 
their results compared. Studies of the possible sensitivities of future measurements at the ILC 
and FCC-ee are presented in Section 5, and Section 6 discusses our conclusions and possible 
directions for future work. 


2 The Covariant Derivative Expansion and the One- 
Loop Effective Action 


The one-loop effective action may be obtained by integrating out directly the heavy particles in 
the path integral using the saddle-point approximation of the functional integral. The contri¬ 
butions to operators involving only light fields can be evaluated by various expansion methods 
for the application of the path integral. Here we follow the Covariant Derivative Expansion 
(CDE), a manifestly gauge-invariant method first introduced in the 1980s by Gaillard [17 and 
Cheyette (18 , and recently applied to the Effective SM (ESM) by Henning, Lu and Murayama 
(HLM) 113] [f The latter piovide, in particular, universal results for operators up to dimension -6 
in the form of a one-loop effective Lagrangian with coefficients evaluated via momentum inte¬ 
grals. This approach applies generally, and greatly simplifies the matching to UV models, since 
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3 For previous analyses, see 
4 We thank Hermes BeluscaTMaito for pointing out to us another recent paper that computes the 


one-loop effective action for certain dimension-6 QCD operators 
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it avoids the necessity of recalculating one-loop Feynman diagrams for every model. However, 
HLM assume a degenerate mass matrix, which may not be the case in general, as for example 
in the ‘natural’ MSSM with light stops. We show here how their results may be extended to 
the non-degenerate case for the one-loop effective Lagrangian terms involved in the dimension -6 
operators affecting the hgg and /177 couplings, with application to the case of non-degenerate 
stops and sbottoms. 

2.1 The Non-Degenerate One-Loop Effective Lagrangian 

We consider a generic Lagrangian consisting of the SM part with complex heavy scalar fields 
arranged in a multiplet 4>, 

Ajv = £sm + (&F(x) + h.c.) + &( P 2 - M 2 - U(x))$ + 0($ 3 ), (2.1) 

where P = iD with D^ the gauge-covariant derivative, F(x) and U(x) are combinations of SM 
fields coupling linearly and quadratically respectively to <3>, and M is a diagonal mass matrix. 
The path integral over <3? may be computed by expanding the action around the minimum with 
respect to <3?, so that the linear terms give the tree-level effective Lagrangian upon substituting 
the equation of motion for <3>: 

£ t f ee = F ] M~ 2 [(P 2 - U)M~ 2 ] n F + 0(4> 3 ), 

n=0 

whereas the quadratic terms are responsible for the one-loop part of the effective Lagrangian. 
After evaluating the functional integral and Fourier transforming to momentum space, this can 
be written in the form 

£ Woo P = i J ^Tr H-(P, - % ) 2 + M 2 + U ) . 

It is convenient, before expanding the logarithm, to shift the momentum using the covariant 
derivative, by inserting factors of e^ Plld l dq,t : 

= ; / (0 Tr Me F « 8 /i> »»(-(P„ - ,„) 2 + M 2 + U)e - p “ a ' e «»]. 

This choice ensures a convergent expansion while the calculation of operators remains manifestly 
gauge-invariant throughout [^] The result is a series involving gauge field strengths, covariant 
derivatives and SM fields encoded in the matrix U(x): 

£floop = i J (04Tr!n [~{G vll d/dq^ + qrf + M 2 + U ], 


where 


n =0 


5 We refer the reader to 
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for technical details and discussions of the CDE method. 
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Here we defined G' Vjl = —iG U)1 with the field strength given by [P u , P^] = —G' v ^. It is convenient 
to group together the terms involving momentum derivatives: 

^lUoop = * J (2^4 ^r + B) j 

where 

A^-U^G^}— + 8U, (2.2) 

B = -q 2 + M 2 + U, 

and we have separated U = U + 5U. 

Expanding the logarithm using the Baker-Campbell-Hausdorff (BCH) formula gives 

ln(A + B) = In (B) + ln(l + H” 1 A) + ^[In B, ln(l + H _1 A)] + ^ [In B, [In B, ln(l + H" 1 A)]] +... 

and, using the identity [In X. Y] = J2n= i where LxY = [X,Y], we see that all 

possible gauge-invariant operators are obtained by evaluating commutators of A and B. 

As an example, we compute the term contributing to the dimension-6 operator affecting 
Higgs production by gluon fusion: 


O g = g 2 \H 2 \G%G a ^. 

The calculation can be organised by writing A as a series in momentum derivatives, 


A = ^A“ 1 - C 


d r 


n =1 


dq ai ...dq Q 


0 d^ 

— _|_ A&i a 2 


dq, 


ai 


dq, 


+ ... 


aiq a2 


2.2 


where each term is obtained by substituting G and U in Eq. 
A " 1 " 2 D —\G' ai ^G' a2 ^ together with the following commutators: 

A 


Here we require only the part 




d 4 q 

( 2 W 


1 


1 


Tr [-B- 2 [B,A] + -B-*[B,[B,A]] . 


We note that M and U are n x n matrices that do not commute in general, which motivates 
the use of the BCH expansion, first applied to the CDE in 18 . Evaluating the commutators 
we find 


^floop * 


d 4 q 

W) A 


1 


Tr \ B- 2 ( ) - -q a q v B 
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3-3 


4 G ^ G 


and using B 1 = —A ^) r) _ 0 (AI7) n , where A = 1 /(q 2 — M 2 ), we see that to obtain operators up 
to dimension 6 requires retaining up to two powers of U, so that we have traces of the form 


Tr(A a UG ,a fl G ,Ufl ) = ]T (A“^G'%Gf , 

i=l 
n n 

Tr(A a U A b U A c G ,a nG'"^) = 

i=i j =i 
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Here we assume G' = diag^^, G' n ) and A = diag(Ai,A n ), where Aj = 1 /(q 2 — m 2 ), and 
U is a general n x n matrix. To evaluate the momentum integrals of arbitrary powers of mixed 
propagators we need to combine them using Feynman parameters: 


d A q 

w? 


q l A“A^ = 


(a + b+ 1 )! 

(a - !)!(£>— 1)! Jo 


r*l 


dZjdZn 


z^z)- 1 


d A q 


( 2 vr; 


l a a+b 

ij 




5(1 — Zi — Zj) 


where A ^ = 1 /(q 2 — m 2 Zi — m 2 Zj ). Taking care in applying the 5 —function in the summation 
over the matrix indices, we finally obtain the following expression valid in the case of a non¬ 
degenerate mass matrix: 


r efr -> 

i-loop (4^)2 




Uij Uji 

2 2 
mfrrij 
1 3 


G' 





(2.3) 


We have checked this result by extending the log-expansion method of |[l3j to the non¬ 
degenerate case by introducing an auxiliary parameter £ and then differentiating under the 
integral sign: 


/^eff 

^1-loop 


d A q 


Tr In [~(G Ufl d/dq^ + qrf + £Af 2 + U] 



where A^ = 1 /(q 2 — £M 2 ) and £ is set to 1 at the end of the calculation. The expansion then 
reads 


reft _ 

^ 1-loop t 


d A q 

( 2 vr? 


d£Tr 



+ U) 


A «M 2 


and yields the same result as in (2.3), demonstrating the consistency of our approach. 

In general the field strength matrix G^ v may not be diagonal, as for example when the 
$ multiplet contains an SU(2)l doublet and singlet, so that we have a 2 x 2 non-diagonal 
sub-matrix VF“ r a involving the weak gauge bosons IF“. The relevant non-degenerate one-loop 

^ ^ * - [TTI 

effective Lagrangian terms then generalise to the universal expression that can be found in 23 u 


2.2 A Light Stop in the hgg and ^77 Couplings 

The result of the CDE expansion is universal in the sense that all the UV information is en¬ 
capsulated in the U, M matrices and the covariant derivative, while the operator coefficients 
are determined by integrals over momenta that are performed once and for all. The simplic¬ 
ity of this approach is illustrated by integrating out stops in the MSSM, whose leading-order 
contribution necessarily appears at one-loop due to R-parity. Since gluon fusion in the SM also 
occurs at one-loop and currently provides the strongest constraint on any dimension-6 operator 
in the Higgs sector, we first calculate its Wilson coefficient within the EFT framework. Later 

6 The previous version of this paper contained an erroneous expression that was correct under the 
assumption of G'^ v commuting with M and U, but not for the fully general case. However this does not 
affect any of our results. 
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we extend the calculation to the the dimension -6 operators contributing to the /177 coupling, 
and comment on the extension to other dimension -6 operators. 

The M and U matrices are given by the quadratic stop term in the MSSM Lagrangian, 


-Cmssm D d> t (M 2 + U(x ))$>, 


where 4> = (Q , t* R ), and 




U = 


(. h 2 t + + y&pHH* - \{gi Y QC 2 p+\<&)w 

htX t H t 


{hi 


h t X t H 

\ 9 lY- tR c 2 p)W 


Here we have defined H = ia 2 H*, lit = VtS/ 3 , X t = A t — //cot/3, and the hypercharges are 
Yq = 1/6,1/^ = —2/3. The mass matrix entries rriQ and m^ R are the soft supersymmetry¬ 
breaking masses in the MSSM Lagrangian. We note that Q = (t^ , 6 /,) is an SU{2)l doublet, 
so U is implicitly a 3 x 3 matrix, and there will be an additional trace over color. Substituting 
this into the CDE expansion with G Rl/ the gluon field strength, we extract from the universal 
one-loop effective action the term 


eff 


^ 1 - 


1o °p D ( 47r )2 24 


1 1 


h l - blew + hi + \glc 2 p 
ml 

tR 


m Q 


This yields the dimension -6 operator O g in the ESM: 


hlX 2 

rn 2 ~m'j 
Q t R , 


g 2 \H\ 2 G^G a ^, 


^dim-6 g Og : 


m 


W 


with the Wilson coefficient given in this normalisation [[] by 


m 


C-a — 


W 


(4?r ) 2 24 


hi ~ \glc 2 p hi + \g\c 2 $ 

O' o 


m 


Q 


m 


tR 


hlX 2 

m 2 -ml 
Q tR , 


This example demonstrates the relative ease with which one may obtain a Wilson coefficient 
at the one-loop level without having to compute Feynman diagrams in both the UV model and 
the EFT that then have to be matched, a process that must be redone every time one adds a 
new particle to integrate out. Here we may add a right-ha nded sbottom simply by enlarging 
the U matrix for = (Q , t* R , b* R ) and plugging it back into (2.3), giving the result 


= _ m w 


1 


(4?r ) 2 24 


h b + h l ~ \ glc 2 p hi + \gjc2p h 2 - \glc 2 p 


h 2 t X 2 


h\Xl 


m 


Q 2 


771.-2 


m b 2 R 


rriAi rrip 


l Q 2 


m 


Q 2 m b 2 


( 2 - 4) 

We compute similarly the dimension -6 operators affecting the /177 coupling, with the field 
strength matrix given in this case by 


G^ = 


W’%T a + Yq 0 iw \ 


-Y t B' 

t-R [IV 


'In general, barred coefficients are related to unbarred ones by c = c^j- where M = v, mw depending 
on the operator normalisation in the Lagrangian. 


7 





















Evaluating this in the CDE method then yields directly 


C, 


. _ cbb rn . cww , CWB 

-dim-6 7) —6 ~ U BB H-n— Vww H- o~OwB , 

1T17,, rn..., m,,.. 


where 
Obb 
and 
(47T ) 2 _ 

-2 ~ C BB 


'W 


m 


w 


w 


m 


'W 


-g 2 \H\ 2 B gu B^ , O ww = g 2 \H\ 2 W« v W a ^ 

1 / 6/q 2 - g\c 2 p 32 [g\c 2 g + 3 hf) \ 

864 \ m\ m~ 


Owe = 2 9l g 2 H^T a HW; u B^ , 


J Q 


tR 


+ h 2 t X 2 


( 


\ 


—103 m%m~ 
Q tR 


— 39 m%m~ + + 16m^ + m~ 

_ Q t R _ Q t R _ Q _ER 

144 m 2 ~rrij ^™ 2 

O t ,1 


( m/^mf 


4m%m~ 1 In / 


r, s"‘i» [ m i - m V ‘ 


+ 


r~rri~ — ^rri^rri- 

Q t R Q t R 




( ~ 771,2 
\ W 


7r/ 


/ 


47t) 2 _ 

9 c ww = 
m w 


6 hi - g\c 2 g 

96 m 2 - 
w 


( 


+ h 2 t X 2 


• + to- 

\ Q tR/ 


V 


(2.5) 

4 ) + m o +” 4 ) ln (^) ^ 

16toT (mT — to- ) 4 4 ( — to- ) 5 

Q V Q t Rj \ Q t Rj y 

( 2 . 6 ) 




CWB = - 


g 2 c 2 p + 2/l 2 

48,4, + 

w 


( 


33— 3m 2 ™ 4 


iX'rrtz ~ orri-^m 
Q t R _ 

24m 2 ~ ( 


j r\t Hr + 5 rrA + 

Q t R Q t R 


m 2 ™ 2 ( 2 m 2 


+ m? ) In 


In the basis used in 1 6 ], the operators Oww and OwB are eliminated and constraints are placed 
on = Obb- The coefficients are related by c 7 = cbb + cww ~ cwb [^J 

To summarise, one may calculate c g and c 7 from integrating out a heavy complex scalar 
<f> in an arbitrary UV model by substituting the SM field matrix, U(x), and field strength 
matrix, G gu , into the CDE expansion. The computation of one-loop Wilson coefficients is thus 
reduced to evaluating the trace of a few matrices. These universal results are extendable to 
all dinrension -6 operators and apply also when integrating out heavy fermions and massive or 
massless gauge bosons [13,23 . 


3 Feynman Diagram Calculations and Comparison 

To estimate quantitatively the validity of the dimension -6 EFT we compare the coefficients 
obtained above with results from an exact one-loop calculation in the MSSM. This is achieved by 
calculating the Feynman diagrams in Fig. [l]then matching the h —» gg and h —> 77 amplitudes 
in the EFT with the equivalent MSSM amplitude. In the EFT the operators O g and 0 7 can be 
expanded after electroweak symmetry breaking (EWSB) around the vacuum expectation value 
v ~ 174 GeV in order to get the Lagrangian 

ChYV = glV2v^hG%G a ^ + gW^v^hB^B^ , 

m w m w 


The log terms in Eqs. [2.5-2.7] cancel in the conribution to c 7 . 



























AOSOS' ah 

h S 1 

— -► 

\ i 

Figure 1: Leading order tree-level Feynman diagram for the EFT (left) and one-loop 
diagrams for the squark contributions (middle and right) to the h —> gg /77 amplitude. 




corresponding to the following Feynman rules for the hgg and /177 vertices: 

h 

m w 


iV hgg(P 2 ’Pk = ~ / Ligl'f 7 iv^ r {p2Pig >lV ~ P2P3) , 
iV/f^(p 2 ,P 3 ) = -^ie 2 V2v-^ r (p 2 p 3 g ftu - P 2 P 3 ) 


w 


Thus the h —>■ gg and h —> 77 amplitudes for on-shell external particles are 


A 

*/v 


hgg _ 


EFT 


= —16g 2 V2v y~ (tUtMl - 2 (£. Pl )(£. Pl )) , 
m w 


A h'y'y _ 
EFT 


= — 2 gf cos 2 9 w V2v^~ {Q.£m£ - 2 (£.pi)(£-pi)) 
m w 


(3.1) 

(3.2) 


where the are the polarization vectors of the gauge bosons. 

We computed the one-loop diagrams in Fig. [l]in the MSSM and checked our results using 
the FeynArts package {27]. The CP-even Higgs bosons are rotated to their physical basis by a 
mixing angle a which we set to be a = /3 — 7 t /2 corresponding to the decoupling limit when the 
pseudo-scalar Higgs mass is much heavier than the mass of the Z gauge boson, as indicated by 
the experimental data 25 and appropriate to our scenario of light stops [^] 

When comparing the EFT and MSSM amplitudes we may choose the momenta of the 
external particles to be on-shell for convenience. The result of this procedure for the h —>• gg 
amplitude yields the same expression as (3.1) with the replacement c g —> c^ SSM , where 


c — = (c R 
^9 '•a 

where the part due to stops is given by 


) + (c“T, 


(3.3) 


77 ? 2 N* 

m W 9 


/—iVloolVlNc _ 

^ 9 > (UA^r\2 


6(4vr ) 2 


a ft = C2/35i 

9 s 

b W 


- v 2 c 2 ggl (2 c 2W + 1 ) + 3 ( 3 v 2 h 2 + 2 ( 'm 2 R - m|) 


C2W + 2 + m 2 R 


+ 36 h 2 [y 2 h 2 t + m| + m\ R - X 2 ) , 


D l = 
9 


v 2 c2ggl 


5 w 


v 2 C 2 fsgl (2 c2w + 1) + 3 ^3 v 2 h 2 + 4 (m 2 R - m 2 kj c 2 w + 4mj~ + 2 m~ R 


+ 36 [v 2 h 2 + 2 rn 2 kj ^ v 2 h 2 + 2 m 2 ^j — 72v 2 h 


2u2 t x?, 


9 The case of relatively heavy stops has been demonstrated to be described in a very compact and 
convenient way, depending only on the two parameters tan (3 and the pseudo-scalar Higgs mass, when 


the observed Higgs mass is taken into account 25 . 
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and the sbottom contribution reads 




m 2 v c 2p9i |6 

[ ml - rnf ) c 2 w + + 2m~ 

LV b R QJ Q b R 

- v 2 c 2 pgl ( c 2 w + 

2)} 

®^ 7r ) 2 ^12m~ — v 1 C2pg\^ 

v 2 c 2 ggl ( c 2 w + 2) - 24m±s( v 




For c 7 we simply have 


-MSSM 

1 


_ _ ^MSSMy 


-MSSM'jb 

2 


+ ) 


(3.4) 


In the limit r->0we obtain the same expressions as c g and c 7 in (2.4) and (2.7), respectively. 
Since c g and c 7 correspond to a truncation of the full theory at the dimension -6 level, they 
contain only the leading-order terms in an expansion in inverse powers of the stop mass, whereas 
the MSSM result is exact and include higher-order terms in v/rrir ? that would be generated 
by higher-dimensional operators in the EFT approach. Therefore, we expect the discrepancy 
between the two approaches to scale with the ratio v/m^~ b for and the differences between 
the EFT and exact MSSM results gives insight into the potential importance of such higher¬ 
dimensional operators. We note that a large value of Xf in terms like v 2 mlyX 2 /m~ could 
potentially affect the validity of the EFT even for large stop masses, but the positivity of the 
lightest physical mass eigenvalue imposes an upper limit Xt — m~/mt- 



Figure 2: Values of AR, defined in (3.11), in the degenerate case jrq = = m, 

for tan/3 = 20 and the indicated values of X t , as a function of mi (left panel), and as 
functions of (right panel). 


The physical mass eigenstates are obtained by diagonalizing the squark mass matrices j24j 


u2 = ( ™ 2 q + ™ 2 LL rn q X q \ 

q V m i x q m \ + m rr ) 

with the various entries defined by 

m LL = m2 Q + ( Iq L — Qq s w ) C2/3 , 

m RR = m \ R + Qq s W M Z c 2 0 i 
X q = A q — g(tan /3)~ 2I ^ L . 


(3.5) 


(3.6) 

(3.7) 

(3.8) 


Q q and l) L is the electromagnetic charge and the weak doublet isospin respectively. After 
rotating the 2 x 2 matrices by an angle 6 q , which transforms the interaction eigenstates (fi and 
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qji into the mass eigenstates q± and < 72 , the mixing angle and physical squark masses are given 
by 


s 2 e q = 


2m q X q 


m 


c 2 e q = 


LL 


— m 


RR 


rn l ~ rn l 


2 2 , x 

m q t ,2 = m q + g 


m LL + m%R =F V {m 2 LL - m 2 RR ) 2 + 4 m 2 X 2 


(3.9) 

(3.10) 


We see that in the stop sector the mixing is strong for large values of the parameter Xf = 
At — pcot /3, which generates a large mass splitting between the two physical mass eigenstates 
and makes q\ much lighter than the other sparticle < 72 - 

We now compare the values of the c g coefficients calculated in the MSSM and the EFT [^J 


A R = 


«EFT 

^n 


^MSSM 


(3.11) 


Fig. |2j displays values of A R for the degenerate case m.Q = m^ R = rrhf, three different values of 
Xt and the representative choice tan/3 = 20. In the left panel we plot A R as functions of m^, 
and the right panel shows A R as functions of the lighter stop mass, . We see that in both 
cases A R <0.1 for rnj(rn^) > 500 GeV, with a couple of exceptions. One is for the relatively 
large value Xt = 3in the left panel, for which A R >0.1 for < 1000 GeV, and the other is 
for Xt = 2777 .^ and rn^ ~ 290 GeV in the right panel, which is due to a node in c^ SSM . These 
results serve as a warning that, although the EFT approach is in general quite reliable for stop 
mass parameters > 500 GeV, care should always be exercised for masses < 1000 GeV. 
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Figure 3: Contours of the differences \(^ ¥r — c^ SSM \ (left panel) and \c^ FT — c^ SSM \ (right 
panel) in (X t /mf,7ni) planes for the degenerate case rriQ = mi R = mi with tan/3 = 20. 
Also shown are contours of mi = 200 GeV, 500 GeV and 1 TeV and regions where the 
t\ becomes tachyonic. 


A similar message is conveyed by Fig. [3j which uses colour-coding to display values of 
the differences |c^ FT — c^ SSM | (left panel) and |c FFT — c^ SSM | (right panel) in ( X t /m^, mf) 

10 We omit RGE effects that mix the coefficients in the running [28| , as they would be higher-order 
corrections beyond the one-loop level of our analysis. 
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planes for the degenerate case nig = m^ R = mf with tan /3 = 20. Also shown are contours of 
rrif = 200 GeV (red), 500 GeV (green) and 1 TeV (yellow) and regions where the t\ becomes 
tachyonic (shaded grey). We see that the differences are generally < 2.5 xlCF 6 for |c^ FT —c^ SSM | 
and < ICE 5 for |c FFT — c^ SSM | when > 500 GeV, even for large values of Xf, but that much 
larger differences are possible for rn^ < 200 GeV, even for small values of Xf. 


4 Constraints on Light Stops from a Global Fit 





Figure 4: Results based on the global fit in ^ , varying c g and c 7 simultaneously but 
setting to zero the coefficients of the other dimension-6 operators contributing to the 
Higgs sector. The dotted, dashed and solid contours on the left denote the allowed 68%, 
95% and 99% CL regions respectively. The middle and right figures show the marginalized 
X 2 functions for c 7 and c g respectively. 


We now discuss the constraints on the lighter stop mass that are imposed by the current 
experimental constraints on the coefficients c g and c 7 , comparing them with the constraints 
imposed by electroweak precision observables via the oblique parameters S and T 1 31 1 , as well 
as the ranges favoured by measurements of the Higgs mass Mh and direct searches at the LHC. 
We note that the S and T parameters are related to the dimension-6 operator coefficients cw, 
cb and ct, as defined in the basis of |6| [j^J through 

S = —- ^( cw + cb) ~ 119 (cw + cb) , 

at(mz) 

T = —-- -ct ~ 129cj'. 

a[m z ) 


We shall quote the electroweak precision constraints on cw + cb and ct instead of S and T, 
in keeping with the EFT approach. The stop contributions to these coefficients were given 
in 112,13 , and Table [I] displays the current experimental constraints on c 9 ,c 7 ,ct and cw + cb 


that we apply. 

The constraints on the coefficients in the penultimate column of Table [l] are taken from a 
recent global analysis 6| of LEP, LHC and Tevatron data on Higgs production and triple-gauge 
couplings. For c g and c 7 we list the current 95% CL ranges after marginalising a two-parameter 


Hn other bases cw and cb may be eliminated in favour of cwb- 
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Coeff. 

Experimental constraints 

95 % CL limit 

deg. m h , 
X t = 0 

Cg 

LHC 

marginalized 

individual 

[-4.5, 2.2] x 10" 5 
[-3.0, 2.5] x 10" 5 

~ 410 GeV 
~ 390 GeV 

c 7 

LHC 

marginalized 

individual 

[-6.5, 2.7] x 10” 4 
[—4.0,2.3] x 10” 4 

~ 215 GeV 
~ 230 GeV 

c T 

LEP 

marginalized 

individual 

[-10,10] x 10~ 4 

[-5,5] x nr 4 

~ 290 GeV 
~ 380 GeV 

Cw + Cb 

LEP 

marginalized 

individual 

[-7,7] x 10” 4 
[-5,5] x 10” 4 

~ 185 GeV 
~ 195 GeV 


Table 1: List of the experimental 95% CL bounds on coefficients used in setting current 
limits on stops, which are taken from |dj/. The marginalized LHC limits are for a two- 
parameter fit allowing c g and c 7 to vary, and the marginalized LEP limits are for a 
two-parameter fit of c T and cw + cb- The corresponding lightest stop mass limits shown 
are for degenerate soft-supersymmetry breaking masses tUq = m^ R = mi with X t = 0. 


fit in which both c g and c 7 are allowed to vary [^J as well as considering the more restrictive 
ranges found when only c g or c 7 % 0 individually, with the other operator coefficients set to 
zero. Similar marginalized and individual 95% CL limits on ct and cw + cb are displayed, 
where the two-parameter fit varying ct and cw + cb simultaneously is equivalent to the S, T 
ellipse, as reproduced in | 6 |. We note that the stop contributions to the coefficients of the other 
relevant operators are far smaller than the ranges of these coefficients that were found in the 
global fit. This indicates that one is justified in setting these other operator coefficients to 
zero when considering bounds on the stop sector, if one assumes that there are no important 
contributions from other possible new physics. 


4.1 Degenerate Stop Masses 

Fig. [5] displays the current constraints in the case of degenerate soft masses mg = mi R = mi 
with aecoupled sbottoms, in the upper panels for as functions of Xt/mi and in the lower 
panels for my 2 as functions of , in both cases for tan/3 = 20. The left panels show the 
stop constraints from the current marginalized 95% bounds on c g (red lines) and c 7 (blue 
lines), and the right panels show the corresponding bounds from the current marginalized 95% 
bounds. The solid (dashed) lines are obtained from an exact one-loop MSSM analysis and 
the EFT approach, respectively. The purple lines show the individual bound from ct in the 
EFT approach. The bounds from cw + cb corresponding to the S parameter are negligible 
and omitted here. The grey shaded regions are excluded because the lighter stop becomes 
tachyonic, and the green shaded regions correspond to 122 GeV< Mh <128 GeV, as calculated 
using FeynHiggs 2.10.3 |29], allowing for a theoretical uncertainty of ±3 GeV and assuming 
that there are no other important MSSM contributions to Mf l . 


12 In any specific model there may be model-dependent correlations between operator coefficients. In 
the case with only light stops and nothing else one expects the relation between c g and c 7 shown in (3.4) 


to hold, as studied in 26 . Here we use the more conservative marginalized ranges shown in the middle 


and right panels of Fig. |4| thereby allowing for additional loop contributions to c g or c 7 . 
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Figure 5: Compilation of the constraints in (upper panels) the (Xt/mf ,mf) plane and 
(lower panels) the (mpmif) plane from (left panels) the marginalized bounds on c g (red 
lines) andc 7 (blue lines), and from (right panels) the individual bounds onc g andc 7 . Also 
shown are the EFT bounds on ct (purple lines), the constraint that the lighter stop should 
not be tachyonic (grey shading) and the region where Mh E (122,128) GeV according to 
a FeynHiggs 2.10.3 calculation assuming no other significant contributions from 
outside the stop sector (green shading). 


We see in the upper panels of Fig. [5] that the c g constraints on rn^ are generally the 
strongest, except for large \X t /mf. We also observe that the MSSM and EFT evaluations 
give rather similar bounds on for \Xf/rnf < 1 and > 2. However, there are significant 
differences for 1 < \XjJrnf < 2, due to the fact that the two evaluations have zeroes at different 
values of Xt/mp The next most sensitive constraints are those from T, parametrised here by 
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Figure 6: Compilation of the constraints in the case of non-degenerate soft mass pa¬ 
rameters, including also sbottom squarks and assuming mg = m^ under the hypotheses 
tan/3 = 20 and X t = 0 (upper left panel), X t = 1 TeV (upper right panel), X t = 3 TeV 
(lower left panel) and X t = (lower right panel). The red (blue) lines show the 

current individual 95% CL constraints from c g (c 7 ) as evaluated exactly in the MSSM 
(solid lines) and in the EFT approach. Additionally, the region compatible with c g is 
shaded pink, the band compatible with Mh is shaded green, and regions disallowed by the 
mixing hypothesis or the appearance of a tachyonic stop are shaded grey. 


the coefficient ct, which become competitive with the c g constraints at large \Xt/m^\, but are 
significantly weaker for small values of Xt/rrif. The constraints from c 7 are weaker still for all 
values of X t /m^, as might have been expected because the global fit in 16] gave constraints on 
c 7 that are weaker than those on c g . Indeed, the c 7 constraint is not significantly stronger than 
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the constraint that the t\ not be tachyonic, as shown by the grey shading in the upper panels 
of Fig. [5j We also note that the LHC measurement of Mh favours \Xf/rrif\ > 2 and values of 
mi that are consistent with the EFT bounds. 

These results are reflected in the lower panels of Fig. [5j where we present the (m^, ) 

planes with the marginalized constraints (left panel) and the individual constraints (right panel). 
The MSSM and EFT implementations of the c g constraint give qualitatively similar results, 
and (except for extreme values of are generally stronger than the constraints from ct, 

which are in turn stronger than the c 7 constraint. We also note that the LHC measurement of 
Mh favours moderate values of m^/m^ and values of rrif or m 1 > 520 GeV. 

The limits on the lightest stop mass for degenerate soft-supersymmetry breaking masses 
m.Q = mi R = Trif with Xt = 0 are shown in the last column of Table jlj 

4.2 Non-Degenerate Stop Masses 

We consider now cases with non-degenerate stop soft mass parameters, allowing also for the 
possibility that the lighter sbottom squark plays a role. We show in Fig. [ 6 ] various planes under 
the hypotheses = rn^ and tan/3 = 20, considering several possibilities for Xf. In all panels, 
the constraints from the individual 95% bound on c g are indicated by red lines and those from 
c 7 are indicated by blue lines (solid for the exact MSSM evaluation and dashed for the EFT 
approach), and the region allowed by the exact calculation is shaded pink. 

The upper left panel is for Xt = 0: we see that in the limit m 1 3> mf the c g constraint 
imposes mi > 300 GeV, with a difference of ~ 20 GeV between the exact and EFT calculations. 
On the other hand, if ru.j 2 = we find raj > 380 GeV, again with the EFT calculation giving 
a bound ~ 20 GeV stronger than the exact MSSM calculation. The corresponding bounds from 
the individual 95% constraint on c 7 are ~ 100 GeV weaker. However, we note that the LHC 
constraint on Mh is not respected anywhere in this plane. 

Turning now to the case Xf = l TeV shown in the upper right panel of Fig. [ 6 j we see a 
grey shaded band around the = mj 2 line that is disallowed by t\ — I 2 mixing, and other 
grey shaded regions where rrij t <C m^ 0 (or vice versa) and the lighter stop is tachyonic. In this 
case the Mh constraint (green shaded band) can be satisfied, with small strips of the parameter 
space ruled out by the c g constraint. The c 7 constraint is unimportant in this case. 

When Xt is increased to 3 TeV, as shown in the lower left panel of Fig. [ 6 j the diagonal band 
forbidden by mixing expands considerably, and the c 7 constraint disappears. In this case the c g 
constraint would allow (m^, rrij 2 ) > (400,1100) GeV on the boundary of the band forbidden by 
the mixing hypothesis, but the Mh constraint is stronger, enforcing (mi , ) > (800,1300) GeV 

along this boundary. 

Finally, we consider in the lower right panel of Fig. [ 6 ] the so-called maximal-mixing hy¬ 
pothesis X t = y/ 6 mi 2 . In this case, almost the entire (rr/ 7] , rrij 2 ) plane is allowed by the 
c g constraint, whereas a triangular region at small rn^ and/or is forbidden by the Mh 
constraint. 

It is interesting to compare the limits on that we find with those found in a recent 
global fit to the pMSSM 130] in which universal third-generation squark masses were assumed 
at the renormalisation scale m ^, the first- and second-generation squark masses were 

assumed to be equal, but allowed to differ from the third-generation mass as were the slep- 
ton masses, arbitrary non-universal gaugino masses ALi. 2,3 were allowed, and the trilinear soft 
supersymmetry-breaking parameter A was assumed to be universal but otherwise free. That 
analysis included LHC, dark matter and flavour constraints, as well as electroweak precision 
observables and Higgs measurements, and found > 400 GeV. The analysis of this paper 
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uses somewhat different assumptions and hence is not directly comparable, but it is interesting 
that the one-loop sensitivity of c g to the stop mass parameters is quite comparable. 


5 Sensitivities of Possible Future Precision Measure¬ 
ments 


Coeff. 

Experimental constraints 

95 % CL limit 

deg. 

x f = o 

mi 

Xt = mil 2 

Cg 

tt /''tll50fb~ 1 
iijV - y 250GeV 

marginalized 

individual 

[-7.7, 7.7] x 10“ 6 
[—7.5,7.5] x 1(T 6 

~ 675 GeV 
~ 680 GeV 

~ 520 GeV 
~ 545 GeV 

FCC-ee 

marginalized 

individual 

[-3.0, 3.0] x 10"“ 
[-3.0, 3.0] x 10 " 6 

~ 1065 GeV 
~ 1065 GeV 

~ 920 GeV 
~ 915 GeV 

c 7 

TT /'rll50fb —1 
iJ ^250GeV 

marginalized 

individual 

[—3.4,3.4] x 10 " 4 
[—3.3,3.3] x 10 " 4 

~ 200 GeV 
~ 200 GeV 

~ 40 GeV 
~ 35 GeV 

FCC-ee 

marginalized 

individual 

[—6.4,6.4] x 10 " 5 
[—6.3,6.3] x 10 " 5 

~ 385 GeV 
~ 390 GeV 

~ 250 GeV 
~ 260 GeV 

c T 

tt /'rll50fb -1 
iJ ^250GeV 

marginalized 

individual 

[- 3 , 3 ] x nr 4 
[—7,7] x 10~ 5 

~ 480 GeV 
~ 930 GeV 

~ 285 GeV 
~ 780 GeV 

FCC-ee 

marginalized 

individual 

[-3,3] x 10 " 5 
[-0.9, 0.9] x 10 " 5 

~ 1410 GeV 
~ 2555 GeV 

~ 1285 GeV 
~ 2460 GeV 

Cw + Cb 

TT /'rll50fb —1 
iJ ^250GeV 

marginalized 

individual 

[- 2 , 2 ] x 10 -4 
[- 6 , 6 ] x 10~ 5 

~ 230 GeV 
~ 340 GeV 

~ 170 GeV 
~ 470 GeV 

FCC-ee 

marginalized 

individual 

[- 2 , 2 ] x 10~ 5 
[- 0 . 8 , 0 . 8 ] x 10" 5 

~ 545 GeV 
~ 830 GeV 

~ 960 GeV 
~ 1590 GeV 


Table 2: List of the 95% CL bounds on EFT operator coefficients from projected con¬ 
straints on Higgs couplings and electroweak precision observables at the future e + e~ col¬ 
liders ILC and FCC-ee. The marginalized limits on c g or c 7 (dr or c-w + cb) are for a 
two-parameter fit allowing c g and c 7 (ct and cw + cb) to vary simultaneously but setting 
other operator coefficients to zero. The corresponding lightest stop mass limits shown are 
for degenerate soft-supersymmetry breaking masses uiq = mi = mi with X t = 0 and 
X t /m- t = 2. 

We saw in the previous Section that the precision of current measurements does not exclude 
in a model-independent way most of the parameter space for a stop below the TeV scale, and 
barely reaches into the region required for a 125 GeV Higgs mass in the MSSM. However, future 
colliders will increase significantly the precision of electroweak and Higgs measurements to the 
level required to challenge seriously the naturalness paradigm and test the MSSM calculations 
of M h . 

In this Section we assess the potential improvements for constraints on a light stop possible 
with future e + e _ colliders. As previously, we perform an analysis in the EFT framework via 
the corresponding bounds on the relevant dimension-6 coefficients, and compare it with the 
exact one-loop MSSM calculation. As representative examples of future e + e~ colliders, we 
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Figure 7: The (X t /mpmf) planes, analogous to those in the upper panels of Fig. [^J 
showing prospective marginalized bounds (left panels) and individual bounds (right panels) 
from the ILC 33l with 1150 fb~ 1 of luminosity at 250 GeV (upper panels) and from FCC- 
ee (32l with 10 4 fb~ 1 of luminosity at 240 GeV (lower panels). In the latter case, the solid 
purple lines are the 95% CL contours for electroweak precision measurements from FCC- 
ee incorporating the projected statistical and systematic experimental errors alone, and 
the dashed purple lines also include theory errors from (351. 


focus on the ILC [33] and FCC-ee 132] (formerly known as TLEP) proposals. The scenarios 
considered here for the ILC and FCC-ee postulate centre-of-mass energies of 250 and 240 GeV 
with luminosities of 1150 fb _1 and 10000 fb 1 , respectively. 

Table [ 2 ] lists the prospective 95% CL limits obtained on c 9 ,c 7 ,cr, and cw + cb from a x 2 
analysis, with the marginalized constraints on c g and c 7 obtained in a two-parameter fit to just 
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Figure 8: Compilation of projected ILC 95 % CL bounds from c g (mf) given by red (blue) 
lines in the vs mi 2 plane, analogous to Fig. [6| with m bi = rn^ and tan/3 = 20. Values 
of X t = 0,1, 3, TeV are shown clockwise from top left. The marginalized limits 

are displayed and the individual bounds are very similar. 


these coefficients, and similarly for ct and cw + cb, corresponding to the T and S parame¬ 
ters respectively, as well as the constraints obtained when each operator coefficient is allowed 
individually to be non-zero. The target precisions on experimental errors for the electroweak 
precision observables m\y , Tz,Ri and Ai at the ILC are given in 133], and those at FCC-ee were 
taken from |32j, and include important systematic uncertainties. The errors on the Higgs asso¬ 
ciated production cross-section times branching ratio are from 34 for the ILC and from 22 


for FCC-ee. The numbers quoted in Table [2] neglect theoretical uncertainties, in order to reflect 
the possible performances of the experiments 0 The treatment of the dimension-6 coefficients 


13 


We also show as dashed purple lines in the FCC-ee panels the weaker constraints obtained using the 
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Figure 9: Compilation of projected FCC-ee 95 % CL bounds from c g (c~ f ) given by red 
(blue) lines in the rn^ vs mf 2 plane, analogous to Fig. ^ with rrr b] = and tan (3 = 20. 
Values of X t = 0,1, 3, sj6rnt 1 m t 2 e V is shown clockwise from top left. The marginalized 
limits are displayed and the individual bounds are very similar. 


in the observables follows a procedure similar to that of the global fit performed in |6j, and we 
use the results of [l4] to rescale the constraint from associated Higgs production. 


5.1 Degenerate Stop Masses 


Contours from possible future constraints on c g , c 7 and ct for the case of degenerate soft masses 
m.Q = mj R = m b are plotted in Fig. |7j using again the value tan f3 = 20. The upper panels 
show results for the ILC, the lower panels for FCC-ee, the left panels show the marginalized 


estimates of theoretical uncertainties in 35 , while noting that these have not been studied in detail. 
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constraints and the right panels show the individual constraints. The grey and green shaded 
regions are the same as in Fig. [7j We see that the marginal and individual sensitivities to 


mr 


from c„ and c~, are very similar, whereas the individual sensitivity of ct are much stronger, 


L, g a iiu l 7 

particularly at FCC-ee. 


We see that ILC is indirectly sensitive to ~ 600 GeV, and that 


FCC-ee is indirectly sensitive to stops in the TeV range. The measurement of the ct coefficient 
at FCC-ee has the highest potential reach, though this will be highly dependent on future 
improvements in reducing theory uncertainties (22,35 . 


The limits on the lightest stop mass for degenerate soft-supersymmetry breaking masses 
uiq = = mi with Xf = 0 and Xf /mi = 2 are shown in the two last columns of Table |2j 

5.2 Non-Degenerate Stop Masses 

Moving on to the non-degenerate case, the c g and c 7 95% CL limits for ILC and FCC-ee are 
plotted in the vs plane for various X t values in Fig. [ 8 ] and [9 respectively. The top 
left, top right, and bottom left plots correspond to Xt = 0,1 and 3 TeV respectively, while 
the bottom right plot is for the maximal-mixing hypothesis Xt = ^6m^ mi 2 . We see that the 
ILC sensitivity to c g begins to probe and potentially exclude parts of the green shaded region 
compatible with the measured while FCC-ee would push the sensitivity of c g constraints 
into the TeV scale. In particular, it could eliminate the entire allowed M^ region for Xt = 3 
TeV. 


6 Conclusions and Prospects 


In light of the SM-like Higgs sector and the current lack of direct evidence for additional degrees 
of freedom beyond the SM, the framework of the Effective SM (ESM) is gaining increasing 
attention as a general framework for characterising the indirect effects of possible new physics 
in a model-independent way. The ESM is simply the SM extended in the way it has always been 
regarded: as an effective field theory supplemented by higher-dimensional operators suppressed 
by the scale of new physics. The leading lepton-number-conserving effects are parametrised 
by dimension -6 operators, whose coefficients are determined by matching to a UV model and 
constrained through their effects on experimental observables. In this paper we have illustrated 
all these steps in the EFT approach for light stops in the MSSM. 

In particular, we employed the CDE method to compute the one-loop effective Lagrangian, 
showing how certain results derived previously under the assumption of a degenerate mass ma¬ 
trix can be generalised to the non-degenerate case. The universal one-loop effective Lagrangian 
can then be used without caveats to obtain directly one-loop Wilson coefficients. The advantage 
of this was demonstrated here in the calculation of the c g and c 7 coefficients. One simply takes 
the mass and U matrices from the quadratic term of the heavy field being integrated out, as 


defined in ( 2 . 1 ), and substitutes it with the corresponding field strength matrix into the uni¬ 


versal expression obtained from the CDE expansion 23 to get the desired operators, without 


having to evaluate any loop integrals or match separate calculations in the UV and EFT. 

Since the hgg and /177 couplings are loop-induced in the SM, the c g and c 7 coefficients are 
currently the most sensitive to light stops. The stop contribution to these coefficients is also 
loop-suppressed, thus lowering the EFT cut-off scale, and it is natural to ask at what point the 
EFT breaks down and the effects of higher-dimensional operators are no longer negligible. We 
addressed this question by comparing the EFT coefficients with a full calculation in the MSSM, 
finding that the disagreement is generally < 10% for a lightest stop mass rn^ > 500 GeV, with 
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the exception of a large \X t \ > or accidental cancellations in the Higgs-stop couplings. 

The constraints on c g and c 7 from a global fit to the current LHC and Tevatron data, 
and the constraints on ct and cw + cr from LEP electroweak precision observables, were then 
translated into the corresponding constraints on the stop masses and Xf . The coefficient c g is 
the most sensitive, followed by ct, which is equivalent to the oblique T parameter. In the case 
of degenerate soft masses, this analysis requires > 410 GeV for X t = 0, and rn^ > 200 
GeV if we also apply the Higgs mass constraint. This is competitive with direct searches and 
is complimentary in the sense that it does not depend on how the stop decays. The limits in 
the non-degenerate case are generally weaker than the Higgs mass requirement, though a few 
strips in the parameter space compatible with Mu can still be excluded. 

The sensitivity of future colliders can greatly improve the reach of indirect constraints into 
the region of parameter space compatible with the observed Higgs mass. The most promising 
measurements will be the hgg coupling and the T parameter, with FCC-ee capable of reaching 
a sensitivity to stop masses above 1 TeV. Thus, FCC-ee measurements will be able to challenge 
the naturalness paradigm in a rather model-independent way. 

As LHC Run 2 gets under way, the question how to interpret any new physics or lack thereof 
will be aided by the systematic approach of the ESM. We have demonstrated this for the case 
of light stops in the MSSM, showing how the EFT framework can simplify both the calculation 
of relevant observables and the application of experimental constraints on these observables, 
giving results similar to exact one-loop calculations in the MSSM. 
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